AD-A103  866  WISCONSIN  UN1V-MADISON  MATHEMATICS  RESEARCH  CENTER  F/G  20/4 

THE  INITIAL  BOUNDARY  VALUE  PROBLEM  FOR  THE  EQUATIONS  OF  MOTION  — ETC(U) 
JUL  81  A  MATSUMMURA *  T  NISHIDA  OAAG29-80-C-0041 

UNCLASSIFIED  MRC-TSR-2237  NL 


ADA103866 


— ^  L'  L.. 


MRC  Technical  Summary  Report  #  2237 


THE  INITIAL  BOUNDARY  VALUE  PROBLEM 
FOR  TOE  EQUATIONS  OF  MOTION  OF 
COMPRESSIBLE  VISCOUS  AND 
HEAT-CONDUCTIVE  FLUID 


Akitaka  Matsumura  and  Takaaki  Nishida 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 

July  1981 

(Received  May  1,  1981) 


DT!C 

ELECTE 

I 

SEP  8  1981 


;ll 


vA 


Approved  for  public  roloose 
Distribution  unlimited 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


i 


National  Science  Foundation 
Washington,  DC  20550 

81  9  08  0?8 


'*9r  >  -  ,.*»i  »<•> 


AVlfcJtl  VJ 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


THE  INITIAL  BOUNDARY  VALUE  PROBLEM  FOR  THE  EQUATIONS 
OF  MOTION  OF  COMPRESSIBLE  VISCOUS  AND 
HEAT-CONDUCTIVE  FLUID 

*  ** 

Akitaka  Matsumura  and  Takaaki  Nishida 

Technical  Summary  Report  #2237 
July  1981 
ABSTRACT 

We  prove  the  global  existence  and  uniqueness  of  solutions  to  the 
equations  of  motion  for  compressible,  viscous  and  heat-conductive  Newtonian 
fluid  in  a  bounded  domain,  with  small  initial  data  and  external  force,  and 
boundary  conditions  of  zero  velocity  and  constant  temperature.  We  also  show 
that  the  solution  decays  exponentially  to  a  unique  equilibrium  state.  The 
proof  uses  an  energy  method  similar  to  the  one  used  in  our  previous  results  on 
the  pure  initial  value  problem  plus  some  new  techniques  for  estimates  near  the 
boundary. 
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SIGNIFICANCE  AND  EXPLANATION 


The  motion  of  compressible,  viscous  and  heat-conductive  Newtonian  fluid 
is  described  by  a  system  of  partial  differential  equations  which  is  of 
hyperbolic-parabolic  type  and  highly  nonlinear.  Having  previously  treated  the 
pure  initial  value  problem  for  these  equations  (see,  e.g. ,  MRC  TSR  #1991, 
2194),  we  now  consider  the  initial  boundary  value  problem  for  motion  within  a 
smooth  bounded  container  whose  walls  are  kept  at  constant  temperature,  subject 
to  a  conservative  external  force,  such  as  gravity.  For  this  problem  we  prove 
the  existence  of  a  unique  smooth  solution  for  all  positive  time,  and  we  show 
that  this  solution  must  decay  exponentially  to  a  unique  equilibrium  state. 
Since  the  system  is  quasi linear  with  respect  to  the  unknowns,  density, 
velocity  and  temperature,  we  need  to  assume  that  the  initial  data  are  close  to 
the  equilibrium  state,  and  that  the  external  force  is  sufficiently  small. 
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§1.  Introduction.  Statement  of  Theorem. 

The  motion  of  a  compressible,  viscous  and  heat-conductive  Newtonian  fluid 


is  described  by  five  conservation  laws: 

r  i 

Pt  +  (pu  5X  =  0  ' 

j 


(1.1) 


i  A  j  i  1 

U  +  U  U  +  —  p 

t  Xj  P  xA  P  Xj  xA 


TfPfu1  +  uj  )  +  M*(uk  )«ij) 


+  fJ 


i  -  1,2,3, 


QpQ  . 

6  +  U39  +  - -  U3  =  - -  ((K0  )  +  f)  , 

t  X.  Pc  X,  Pc  X.  X. 

j  V  j  V  j  3 


where  x  -  (x1,x2,x3)  e  R3,  t  >  0,  p  is  the  density,  u  *  (u3,u2,u3)  is  the 

velocity,  0  is  the  absolute  temperature,  p  =  p(P,6)  is  the  pressure,  f  «■ 

(f1,f2,f3)  is  the  external  force,  P  =  P(P,0)  and  P’  =  P'(P,0)  are 

viscosity  coefficients,  <  =  <(P,0)  is  the  coefficient  of  heat  conductivity, 

c„  *  c  (P,0)  is  the  heat  capacity  at  constant  volume  and  'f  *  —  (u3  +  uk 

vv  2  x.  x 

12  31  ^  o 

+  P'(uJ  )  is  the  dissipation  function.  Let  ft  be  a  bounded  open  set  in  R 

Xj 

with  a  smooth  boundary  9ft.  We  consider  the  initial  boundary  value  problem 
for  (1.1)  in  ft  with  the  boundary  conditions 
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(1.2) 


u(t,x)  ■  0,  0{t,x)  »  x  6  30,  t  >  0  , 


where  "9  is  any  fixed  positive  constant;  and  with  initial  data 


(1.3) 


(P,u,0)(O,x)  -  ( Pg #Ug , ) ( x) ,  x  e  0  . 


We  seek  a  classical  solution  for  (1.1)  '  (1.3)  globally  in  time. 

Historically,  a  local  solution  in  time  for  (1.1)  '  (1.3)  was  constructed  by 
Tani  (6].  On  the  other  hand,  Kazhikhov  and  Shelukhin  [3]  proved  the  existence 
and  uniqueness  of  a  global  solution  in  time  for  the  one-space-dimensional 
model  of  (1.1)  '  (1.3).  In  this  paper,  we  prove  the  existence  and  uniqueness 
of  a  global  solution  in  time  for  the  original  problem  (1.1)  '  (1.3)  if  the 
initial  data  and  external  force  are  suitably  small.  Although  the  proof  is 
given  by  an  energy  method  similar  to  those  of  our  previous  papers  [4]  and  [5] 
on  the  initial  value  problem,  the  initial  boundary  value  problem  requires  now 
a  priori  estimates  of  the  solution  near  the  boundary  dft.  In  order  to  state 
the  main  theorem  precisely,  let  us  list  our  assumptions: 

(i)  The  external  force  f*  is  generated  by  a  potential  function  *(x), 
i.e.. 


(ii)  P,  P',  K,  P,  cy  e  C  (O),  where 

O  -  {(P,0)  |  (P,0)  e  (0,+»)  x  (0,+»)). 

( iii)  P,  *,  p,  pp,  Pg,  cy  >  0  and  P’  +  j  P  >  0  for  (P,0)  e  O. 

(iv)  (PQ,  uQ,  0Q)  6  H4(fl),  ♦  e  H5(fl)  and  (Po,0Q)(x)  >  0  for  x  e  ff. 

(v)  (compatibility  condition) 


*  ♦  u_  ^  )  ♦  Plu„  *  6*3)  -  P.4 


- (p_ )  +  (Pn(U-  ♦  un  J  )  ♦  P'u. 

0  x  0  OiX  OsX  0  0#x. 

i  j  i  K 


*  '"oVxj’xj  *  Mail  ■  0  - 


x  Ox 

3  31  an 


-  o  , 


Ml 


where  pQ  -  P^PO*0O^»  u0  *  U*Po'90* ' *  *  * '  and  80  on*  Here  ”  <a)  denotes  the 

Sobolev's  space  on  ft  with  the  norm  1*1^,  end  Hg(ft)  denotes  the 

completion  of  Cg(ft)  in  11^(8)  •  Define  1*1  ^  by 

CK 


and  also  define  P  by 


lul  »  max  sup  | (t“)  u(x)|  , 

C*  lal <k  xeft 


p  "  l  Po(x)dx  ' 


where  V(Q)  represents  the  volume  of  ft.  We  say  that  (p(x),  u(x),  9(x))  e 
1  -  2  -  2 

C  (ft)  *  (C  (ft))  is  an  equilibrium  state  of  the  problem  (1.1)  "v  (1.3)  when 

v  v  v 

(P,u,9)  satisfies  (1.1)  and  the  following  additional  conditions 


(1.4) 


u,3ft  "  0  *  0,3ft 


0  '  ,77n7  /  P  dx 


Then  we  have  the  following  main  theorem. 

Theorem  1.1.  Under  the  assumptions  (i)  "*  (v),  there  exist  positive 
constants  e,  a  and  C  such  that  if  Ip  -P,  u„,  9-01  +  4  e,  the 

U  0  0  4  5 

problem  (1.1)  ”  (1.3)  has  a  unique  global  solution  in  time  (P,u,9)  and  a 

V  V  V  V  __ 

unique  equilibrium  state  (P,u,9)  -  (P,0,9)  satisfying 

p  e  c°(0,+*j  H4(ft))  n  c1(o,+«,»  H3(ft))  , 

(u,9)  e  C°(0,+«;  H4(ft))  n  H2(ft))  , 


«(P,u,9)(t)  -  (P, 0,0)1  <  Ce" 

c 


§2.  Existence  and  Uniqueness  of  the  Equilibrium  State. 

Let  us  rewrite  the  equations  and  conditions  for  the  equilibrium  state 


(P»u, 8) : 


(2.1) 


(PuJ)  »  0  , 

Xj 


(2.2) 


(2.3) 


(2.4) 


’■4  v  v  v  vj  v4  v  n4*  i  4 

Ju  +  p  +  p$  -  ( u(u  +  uJ  )  +  M'u  6  J)  =  0  , 

x.  x.  x.  x,  x.  x,  x. 

j  1  1  j  A  X  3 


V  N*  V  w  V 

Pc  u30  +  0pflr  -  (X0  )  -  ¥  =  0  , 

v  Xj  ^0  xj  X.  X. 


ul  3JJ  =•  0  ,  9|  ao  “  0  ,  /  p  dx  *  /  P  dx  , 


30 


where  p  =  p(P,9),  P  ■  p(p,0),...#  and  so  on.  Then  we  have 

Lemma  2.1.  Given  the  assumptions  (i)  ~  (v),  there  exist  positive  constants  e 

and  C  such  that  if  l®l  4  e,  the  problem  (2.1)  ~  (2.4)  has  a  unique 
_  C  _  _ 

solution  (P(x),O,0)  in  a  small  neighborhood  of  (P,0,9)  in 

0^(0)  x  (C2(0))2  satisfying 

Ip  -  pi  i  <  cl*l  1  , 

C1  C1 

where  P(x)  is  determined  by  ^  _ 

r  PD(s'®> 

Const.  -  *(x)  =  /  - 


(2.5) 


ds  , 


/  p(x)  -  P  dx  *  0 

0 


Furthermore  we  have  P  e  H  (0)  and 

(2.6)  Ip  -  pi  4  Cl*l  for  l*l_  4  e  . 

5  5  5 

Proof  of  Lemma  2.1.  Suppose  |p-p|,  |u| ,  | 0—6 |  <  jmin(P,0).  Then  we 
may  estimate  the  equalities 


r 


P  p  (s,6) 

/  [2.1]  x  /  -E -  ds  dx  =  0  , 


(2.7) 


<  /  (PuA) [2.2] Adx  *  0  , 


/  P(0-0) [2.3]dx  =  0  , 


V 


-4- 


where  [2.1],  [2.2]*  and  [2.3]  denote  the  left  hand  side  of  (2.1),  (2.2)*  and 
(2.3)  respectively.  After  integrating  (2.7)  by  parts  and  using  Poincare's 
inequalities 

(2.8)  lul  <  ClDul  ,  80-91  <  C<D0l  , 

one  can  show  that 

(2.9) 


iDul2,  Id9I2  <  C(lDpi  +  80-91  )(«Dul2  +  «D0l2)  . 

C°  C° 


Here  and  in  what  follows,  we  denote  by  D  f  for  k  >  0  and  z  -  (z.,...,z  ) 

z  i  in 

the  vector  function 

a  a  a 

D*f  *  O^f/dz.1  3z  2  ••••  3z  m  for  all  |a|  =  k) 
z  12  m 


and  in  particular  write  DKf  when  z  =  x.  Therefore,  if  lDpl  and 
_  C° 

•®-08  are  suitably  small,  it  follows  from  (2.4)  and  (2.9)  that 

C  .  _ 

(2.10)  u  =  0  ,  650  . 

Substituting  (2.10)  into  (2.2),  we  have 


(2.11) 


P  p(s,9) 

(/  —  ds  +  *)  =0  * 


which  implies  (2.5)  immediately. 

One  may  use  the  implicit  function  theorem  once  to  see  that  the  constant 

5 

in  (2.5)  depends  on  ®  as  a  smooth  functional  on  H  (H),  and  again  to  prove 

(2.6). 


§3.  Local  and  Global  Existence. 

First  let  us  rewrite  the  problem  (1.1)  ~  (1.3)  by  the  change  of  variables 
(P,u,9)  ♦  (P+P,  u,  9+9)  as  follows: 


(3.1) 


p  +  u3P  +  (p+p)u3  +  P  u3  «  0  , 

t  Xj  Xj  Xj 

u*  +  u3u* - ^3  (Ulu*  +  u3  )  +  U' yk  0i3)  + 

*  Xj  P+P  Xj  Xi  ’Sc  Xj 

pp  p0  Pp  (P+P,e+?> 

♦  -s  ♦  —  •. .  ^  - 1)»  . 


P+P  i  p+p  i  <P+P)pp(P,9) 


X. 

l 


(0+*>P0 

0  +  u30  +  - 3 — -  u3  - 

^  *  Xj  (P+P)c  Xj  (P+P)c 


1 


(Of9v  >*  +  *>  , 
j  J 


(u,0)|3n  -  0  , 


(3.2) 

(3.3) 

Furthermore  we  rewrite  the  problem  (3.1)  ~  (3.3)  as  follows* 


(P,u,9) (x,0)  -  (Po'uo'%^x)  * 


(3.4)' 


L°(p,u)  =  P  +  u3P  +  P  u3  -  f°  , 
u  t  Xj  Xj 


( 3.4) 1  lNp.u,©)  a  u*  -  y  U1  -  (y  +  y')u3 


Vi 


x.  x. 
1  j 


p.p  +  P-0  ■  fx  , 

1  X.  2  X.  ' 


(3.4) 


L4(u,0)  =  0  -  k0  +  p  u3 

t  x^ Xj  ^3  x_ 


(3.5)  <u'e> 1  an  “  °»  <p.u,0)(o)  -  <VvV  r 

where  define  L^(P,u)  =  P.  +  v3p  +  Pu3  and  M  *  y/P,  y'  *  y'> 
v  t  Xj  Xj 

Pi  =  Pp/P#  p2  *  p0/p»  p3  *  ®  pp /p  cy,  <  =*  k/p  cv,  y  -  k(p,0),  u' 

y'(p,0},  Pp  *  Pp(P»0)  and  so  on,  and 

f°(P,u,0)  =  (p-p-p)u3  -  P  u3  , 

XJ  *J 


P  e  c°{t1,t2>  h4)  n  L2(t1,t2»  h4)  , 

Pt  e  c°(t1,t2»  h3)  n  l2 ( 1 1 »t2 1  h3)  , 

Ptt  e  c°(t1ft2;  H1)  n  L2(t1,t2>  h2)  , 

(u,e)  e  c°(t1,t2>  H4  n  hJ)  n  t.2(t1ft  i  h5)  , 

(u,8)t  e  c°(t1,t2»  H2  n  h^)  n  L2(t1#t2»  H3)  , 


-7- 


I 


(u,9)tt  e  c°(t1,t2»  l2)  n  L2(tlft2»  Hq)  , 


N(t1#t2)  <  E>  . 


Here  is  defined  by 

N2(t.,t  )  »  sup  (*P(t)*2  +  *P  (t)*2  +  Ip  (t)l2  + 

1  2  t1<t<t2 

+  l(u,9)(t)l2  +  *(u  , 9  )(t )  I2  +  Mu  ,9  )(t)l2)  + 

4  t  t  Z  tt  tt 


+  I.2  (IP(T)»2  +  IP. (T)l2  +  *P  ( T) I2  +  I(U,9)(T)I2  + 
t1  4  t  3  tt  2  5 


+  *<u.,9  )(T)I2  +  Mu  ,  9  )  ( T)  l2)dT  . 

t  t  3  tt  tt  1 


We  can  get  a  global  solution  by  a  combination  of  a  local  existence  result  and 
some  a  priori  estimates  for  the  solution  in  X. 

Proposition  3.1.  (local  existence)  Suppose  the  problem  (3.1)  '■  (3.3)  has 
a  unique  solution  (P,u,9)  e  X(0,h|  Eg)  for  some  h  >  0  and  then  consider 
the  problem  (3.1)  ~  (3.3)  for  t  >  h.  Then  there  exist  positive  constants 


f ,  6q  and  Cg  (eQ  /  1+Cg  <  Eg)  which  are  independent  of  h  such  that  if 
N(h,h),  191  <  eQ,  the  problem  has  a  unique  solution 

(P,u,9)  e  X(h,h+T;  CgN(h,h))  . 

Proposition  3.2.  (a  priori  estimate)  Suppxsse  the  problems  (3.1)  ~  (3.3) 

has  a  solution  (P,u,9)  e  X(0,h>  EQ)  for  some  h  >  0.  Then  there  exist 

positive  constants  e1  and  (ei  <  eg,  e^^  <  EQ)  which  are  independent 

of  h  such  that  if  N(0,h),  191,.  <  e^,  it  holds  that  N(0,h)  <  €^(0,0). 

Given  Propositions  3. •  and  3.2,  the  global  existence  uniqueness  of 

solutions  may  be  proved  as  follows.  Choose  the  initial  date  (P^,u  ,9  )  and 

0  0  0 

the  potential  function  9  so  small  that 


t 


*  .*»%.  <*•*** mft 


. 


I 


I 


i 


l 


N(  0,0 )  <  min(e  ,  ~r,  - ,  1  . ’5)  * 

0  c/  1  +  co 

and  I#»  <  Then  Proposition  3.1  with  h  =  0  gives  a  local  solution 

(P,u,0)  e  X(0,t,  CqN(0,0))  . 

Since  CqN{ 0,0)  <  e.]  <  eQ,  Proposition  3.2  with  h  =  T  implies 

N(0,T)  <  0^(0, 0)  , 

and  Proposition  3.1  with  h  =  T  implies 

(P,U,9)  e  X(T,2T;  CqN(t,t))  , 

e  X(0,2T;  /  I+Cq  N(0,t))  . 

Hence,  since  /  1+cJ  N(0,T)  <  c^/l+cj  N(0,0)  <  Proposition  3.2  with 

h  =  2T  implies 

N(0,2T)  <  0^(0, 0)  , 

and  Proposition  3.1  with  h  =  2t  gives 

(P ,U, 0 )  e  X(2T,3T;  CQH(2T,2T))  , 

e  X(0,3T,  /  1+Cg  N( 0, 2T) )  . 

Repetition  of  this  process  yields: 

Proposition  3.3  (global  existence) 

There  exist  positive  constants  e  and  C  such  that  if  N(0,0), 
l*l5  <  e,  the  problem  (3.1)  ~  (3.3)  has  a  unique  solution 

<P,u,9)  e  X(0,+“r  CN(0,0) )  . 

§4.  Proof  of  Local  Existence  and  Uniqueness. 

We  note  only  how  to  prove  local  existence  and  uniqueness  for  h  =  0,  and 
omit  the  details.  Suppose  (h,v,C)  e  X(0,t;  eq).  Then  we  first  need  to  solve 
the  following  equations: 


-9- 


(4.1) °  L°(P,u)  -  f°(h,v,C)  , 

(4. 1) 1  -  f1(n/v/o  , 

(4.4)4  L4(u,fl)  -  f4(n,v,C)  , 

(4.2)  <u,e)i3Q  “  °'  cp,u»«)(o)  -  <Vuo'V  ' 

where  we  note  that  (f1)4^  satisfy  the  conditions 

f°  e  C°(0,T;  h3)  r>L2(0,t,  H4)  , 

f°  ec°(o,t,  t2)  ni2(0(i|  H2)  , 

4  e  iW,  i2>  , 

fA  e  c°(o,t,  h2)  m2(0,i|  h3)  , 

f*  e  c°(o,x,  t2)  n  l2<o,t,  h1)  , 

4  6  C°(0'T’  H’1)  '  1  <  1  4  4  ' 

where  H-1  represents  the  dual  space  of  Hg. 

To  solve  (4.1),  (4.2),  we  may  regard  the  equations  (4.1)  as 

(4.2)0  p  +  vjP  -  -P  +  f°  , 

t  Xj  Xj 
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(4.2)- 


i  A  i  A  A  4  i  . 

u.  -  V  u  -  (U+U')uJ  ■  f  -  p_8  -  p„p 

*  xjXj  Vj  *  Xl  P1  Xl 


(4.2) 


8  -  <  8  »  f  *  -  p  u^  , 

t  XjXj  j  Xj 


and  use  standard  arguments  about  first-order  hyperbolic  equations  for  (4.2)° 

and  second-order  parabolic  systems  for  {(4.2)*}*^ .  Once  we  get  the  solution 

of  the  linear  problem  (4. 1 ) (4.2 ) ,  we  may  construct  an  approximation  sequence 

{(P,u,0) )+  .  for  the  nonlinear  problem  (3.4) (3.5)  as  follows: 
n~o 

(p.u,0)(O)  s  (P0,u0,e0)  , 

and  for  n  >  1, 


r 


T0  (n)  (n),  .0.  (n-1)  (n-1)  B(n-1). 

L  (n_i)(p  »u  )  *  *  (P  •  «  *  0  >  # 


i .  (n)  (n)  «(n)  #i ,„<n-1)  (n-1)  B(n-1>, 

h  (P  ,  u  #8  )  •  f  (P  ,  u  ,0  )  , 


_ 4.  (n)  fl(n) .  4. ft(n-1>  ( n-1)  B(n-1). 

l>  (u  ,8  )  *  F  (P  ,  U  ,8  )  , 


v  (u(n),  0<n))l3n-  0  »  <p(n).  u(n),  ®(n))(0)  -  <VvV 


Finally,  we  may  apply  the  contraction  mapping  principle  to  the  nonlinear 
problem  (3.4),  (3.5)  in  the  space  X  under  the  smallness  conditions  on  the 
initial  data  and  external  force,  to  obtain  existence  and  uniqueness  of  the 
solution. 
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§5.  Proof  of  the  A  Priori  Estimates. 

Throughout  this  section,  we  suppose  that  *41^  <  EQ  and  that  the  problem 

(3.1)  '•  (3.3)  has  a  solution  (P,u,8)  e  x(0,h;  E^ )  and  we  note  that  all 
constants  are  independent  of  h.  First,  recall  Poincare's  inequalities  for 

(u,8); 

(5.1)  lul2  <  ClDul2  ,  I8I2  <  ClD6l2  , 


and  the  following  lemma  for  P. 
Lemma  5.1.  It  holds  that 


Ipl2  <  ClDP*2  . 


This  is  proved  using  the  fact  that  the  equation  (3.1)°  implies 

/  P(t)dx  -  /  PQ  dx  =*  0 

for  all  t  >  0  (cf.  [7] ). 

Next,  since  the  viscosity  and  heat  conduction  terms  of  (3.4)  define 

k  1 

strongly  elliptic  operators  on  (u,8),  we  have  for  (u,9)  e  H  n  HQ 
(1  <  k  <  5)  (cf.  11],  [2] ): 


(5.2) 


lul2  <  C  I  l5  +  (P+Q' )u^  I?  _  , 

k  i-1  *jxj  xi  j  k"2 


i  o i2  <  ciic  e  i2  . 

k  c  xjXj  k-2 

Equations  (3.4)  and  estimates  (5.1)  and  (5.2)  easily  yield  the  following 
1 6nDU)  # 

Lemma  5*2#  For  1  <  k  <  5#  it  holds  that 

"'i-  *  K-,  ♦  j,  - 

i-1 

,9,J  ‘  '"'t'w  *  '<i  ♦  • 


Third,  the  following  estimates  on  the  time  derivatives  of  the  solution  over 
all  of  fl  can  be  obtained  rather  easily. 


Lemma  5.3.  For  any  positive  number  e  and  0  <  k  <  2,  it  holds  that 


lD*(P,u,9)(t)lZ  +  J  Id£(u,9)IZ  +  lD^(~)l2dT 


<  dD*(P,u,6)(0)lZ  +  /  e  iD^pl2  +  C(1  +  e_1)(iD*fuiz  + 


.“1 . .  «_k_0,2 


+  l  iD^f1!2  )dT  +  CN3(0,t) 
i=1  t 


where  —  is  defined  by 
at 


dp  3p  j  3p 
4t  *^+  “  ^  * 

This  lemma  is  proved  by  estimating  the  following  equation  derived  from  (3.4): 


r  k.,0  .  _k„  .  i  .  i  ,  P2  k  4  kft 
J  J  —  D.  L  •  DP  +  D.  L  •  D.  u  +  —  D.L  *  D  B  dxdT 
0pfc  t  t  t  p3t  t 

t  P 


f  1  k  0  K  ,  ^X.l  K  1  ,  l  _k_4  .  _Kft  . 

J  —  D  f  •  D.  P  +  D.  f  •  D.u  +  —  D  f  •  D  0  dxdT  , 
'  ~  t  t  t  t  P„  t  t 


k.i  .  k  i  /2  k,4 


0  P 


and  noting  from  (3.4)  that 


„J  ..is,’,*  . 

-j  5*  p 

Fourth,  let  us  establish  the  interior  estimates.  Let  XQ  he  any  fixed 
function  in  CQ(0).  Then  we  have 


Lemma  5.4.  For  0  <  t  4  3,  it  holds  that 


»X0Dt+1P(t)  «2  +  /  Ix0D*+1P*2<3*  < 

0 


M 1-2. 


«  *  vi  *  Mt,,«  *  * 
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+  ,f0,Li +  I  ,fi,5 dT  +  . 

*+1  i-i  * 


Lemma  5.5.  For  any  positive  number  e  and  1  <  k  <  4,  it  holds  that 


IX0Dk(P,u#e){t)l2  ♦  /  lx0Dk+1u»2  ♦  *X0Dk(||)»2dT 


<  c»p0,u0,6oi£  +  /  e  »x0Dkpi2  ♦  c(i  ♦  e"S  ( 


*u,8*2  +  lf°l“  +  l  lfilj_1)dt  ♦  C»3<0,t)  . 


To  prove  these,  we  may  estimate  for  Leona  5.4  the  equation 


L  1 •  -V  *  ~k  *>v  • 


i  2M+JI' 


P  dxdT 
Xi 


I  I  XflD4?0  •  D*P  +  — P-~  X^D^f1  •  D*P  dxdT  , 
0  U  Xi  Xi  2U+111  0  Xi 


and  for  Lemma  5.5,  the  equation 


f*  f  fi  v2nkT°  .  ■  2  k  i  (  k  i  .  P2X0  k  4  .  kft  _ 

J  J  —  XaD  L  •  D  P  +  X-D  L  •  D  U  +  -  D  L  •  D  9  dxdT 

0  P  °  0  P3 

f  f  P1  2  0  k  .  2  k.i  lhkli  P2X0  Jc.4  .  k„  . 

J  J  —  X„Df  *  D  P  ♦  XnD  f  •  D  u  +  -  D  f  •  D  8  dxdT 

0  p  °  0  p3 


Fifth,  let  us  establish  the  estimates  near  the  boundary.  We  choose  a  finite 
number  of  bounded  open  sets  (o. .  in  R3  such  that 


and  in  each  open  set  Oj  we  choose  local  coordinates  (+,v,r)  as  follows: 
(i)  the  surface  0^  n  30  is  the  image  of  smooth  functions 
y*  *  y*  (♦/¥>)  satisfying 


(5.3) 


•  Yj  •  i  »  i 

y*  *  y„  >  4  >  °  » 


Y+  *  y*  -  o 


where  6  is  some  positive  constant  independent  of  1  <  j  <  R, 


(ii)  any  x*  in  0^  are  represented  by 


(5.4) 


x1  -  x*(t,i,r)  -  ra1^,*)  +  yi(t,«)  , 


where  n  (♦,<?)  represents  the  internal  unit  normal  vector  at  the  point  of  the 
surface  coordinated  ( ♦,  <fi) . 

Here  and  in  what  follows,  we  omit  the  suffix  j  for  simplicity.  Let  us 

define  the  unit  vectors  e*  and  e*  by 

i  i  i  1,,  i, 

e1  *  y*  '  *2  “  V(y  1  * 

Then  note  that  there  exist  smooth  functions  of  ($,<?),  (<*#0,Y )(♦,*)  and 
(<*,8*  ,Y’ )  (♦,<?)  satisfying 


0  -Y  -a 


Y  0  -0 


a  0  0 


0  -Y’  -a' 


Y»  o  -8 


a*  0*  o 


t5 


An  elementary  calculation  will  show  that  the  Jacobian  J  of  the 
transformation  (5.4)  is 


(5.6) 


J  “  •**  x  V 

-  Iy*l  +  +  S’)r  +  (<*B>  -  Ba')r2  . 


By  (5.6) 

needed. 

as 


,  we  can  see  the  transformation  (5.4)  is  regular  choosing 

Therefore  the  functions  (»M#r)  (x)  make  sense  and  are 

*i 


(5.7) 


\  ‘  5  "V  *  “t’l  '  J  '“I  *  "i1  ' 


%  -  I  "S:  *  Vl  ’  J  <C*1  *  “s’  ■ 


K  m  j  (x* x  Vi  “  ni  ' 


where  A  -  ly^l  +  B’r,  B  -  -ra*,  c  -  -8r,  D  -  1  +  ar  and  J  -  AD 
Hence  (5.7)  implies 


1  i  i  9  i  i  13  i3 
j  (Ae1  +  Be2)  J*  +  J  (Ce1  +  De2)  ft  +  n  K 


Thus,  in  each  Oj,  we  can  rewrite  the  equations  {(4.1)*}|-Q  in 
coordinates  (<l>,^»r)  as  follows: 


Pt  ♦  -  ((Ae^n1  +  Be^uSp^  +  (Ce*  +  De^u^  +  u*n  JPy) 
♦  J  ((Ae*  +  Be*)^  +  (Ce*  +  De*)u*  +  “  f° 

A 

\  *  h  (<a2  +  b2)uh  +  2(AC  +  BD,v  +  (c2  +  °2)ul* 

j 


* 


r  small  if 
calculated 


-  BC  >  0. 


the  local 
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2  1  » 

+  J  u  J  +  first  order  terms  of 
rr 


u  + 


A  A 

1  i  i  u+u*  dp 

+  a  (Aei  +  +  W 


1  i  i  U+ft'  dp 

♦ j  «*; ♦  £ + w + 


A  A 

t  i,u+Hl  dp  ^ 

*  "  'T"  *  P'C>r 

P 


A  A. 

f1  -  P  9  +H£1«»  , 

2  xi  p  xi 


Where  we  note  that  J2  -  (AC  +  BD)2  -  (A2  +  B2)(C2  +  D2).  Now  let  X  be  any 

00 

fixed  function  in  CQ (O^ ) .  By  methods  similar  to  those  of  Lemma  5.5,  we  can 


Lemma  5.6.  For  any  positive  number  8  and  1  <  k  <  4,  it  holds  that 

,x  Di,v>  <p'u)(t)|2  +  f  **  DvDu,Z  +  **  D*^(^,,2dT 


<  c\'»0'2k  *  /  e  *X  d*^p»2  +  c<*  +  e-1)  < 


«u,8l2  +  If0 12  +  l  If1 12  ^JdT  +  CN3(0,t)  . 


Moreover,  by  (5.8),  we  gets 


(5.9) 


A  A. 

2U+U‘  . dp.  . 

“  dt  r  +  p1Pr 


i  i  0 
-n  u  +  —  ( 

*  J2 


(A2  ♦  B2)niu^  +  2 (AC  +  B0)nltt^  +  (C2  +  D^nSi^, 


-  J(Ae^  +  Be*)u^.  -  J(Ce*  +  De*)u^)  + 
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*  *-  '«•  dura— 


,  MA  A.  M  _  »  ,  ,  i-i  ^  2&+M*  -0 

+  first  order  terms  of  (u,0)  +  n  f  ♦  — — —  f 

P 


Estimating  the  equation 


/  /  X2  |D^  D*(equation  (5.9))|2dxdt  , 


for  0  <  k+i  <  3,  we  have 

Lemma  5.7.  For  0  <  k+*  <  3,  it  holds  that 


*  'x  -vr > >,2<”  ‘  * 


♦ c  i  '*  OX  *  \'lt  *  ‘"-eik*i*i  * 


♦  *  i  ,'i,^tdT  *  o,3<°-t>  • 

i-i 


Sixth,  we  note  the  following  lemma  about  the  linear  stationary  Stokes  Equation 
(see  17],  for  example). 

Lemma  5.8.  Consider  the  problem 


f  P  nl  -  ,°  . 


(5.10) 


-  *  jXj  +  Vx4  -  **  • 
u,a a  m  0  ' 


where  g^  8  and  g^  8  H*  (k  *  0). 


k+1 


Then  the  problem  (5.10)  has  a  solution  (P,u)  such  that  P  8  H  and 

1(42  1 

u  8  H  n  H0  which  is  unique  modulo  a  constant  of  integration  for  Pi  this 


solution  satisfies 


‘  c»*0,lt  *  j,  • 

13  ^  * 

Apply  Lemma  5.8  to  the  equations  {(3.4)  )^-0#  we  have 
Lemma  5.9.  For  2  <  k  <  5,  it  holds  that 


'-'l  *  *»0,*-2  ‘  *  ViU  * 


♦  ,Mk-i  *  J,  •» 'li  ♦ 


Furthermore  by  differentiating  the  equations  {(3.4)A}^0  X  times  with 
respect  to  1*  and  V  in  each  0 j ,  multiplying  them  by  x  and  again 
applying  Lemma  5.8,  we  obtain 


Lemma  5.10.  For  1  <  44k  <  3,  it  holds  that 


.  2+4  k  .2  ,  _1+4  k  .2 

■x  D  D^u*  +  *x  D  D*^pl 


« «'*  ■’,X»,£"2  *  *  •<*  * 


*  '°A*  ♦  ,f0,Ui  *  *  ,tl,h'  ■ 

i-i 


Using  all  of  the  above  preparations,  we  are  ready  to  establish  the 
desired  a  priori  estimates.  Although  we  omit  the  details,  by  combining  Lemma 


5.3  with  k  -  0  and  1,  Lemma  5.2  with  k  -  1,  2  and  3, 


5.4  with 


4-0  and  1,  Lemma  5.5  with  k  -  1  and  2,  Lemma  5.6  with  k  -  1  and  2, 
Lemma  5.7  with  k  -  4  -  0  and  k+4  -  1,  Lemma  5.9  with  k  -  2  and  3  and 
Lemma  5.10  with  4-0  and  k  -  1,  and  by  choosing  6  suitably  small,  we 


(5.ii)  ip,u,ei^  ♦  +  ,Pt*2  + 


19' 


*  r  •».«’  *  >»i  *  >£ >2  *  <vv?  * 


<  Clip  ,u  ,e  r  +  sup  (lf°l2  +  l  If1!2) 

2  0<T<t  i=1 


+  /  »f° *2  +  !  if1',  +  *f°'2  +  !  »f^»21dr) 

0  2  i-1  i-1  z  ~ 


+  CN  (0,t)  . 

In  the  same  way,  after  considering  the  situation  differentiated  with  respect 
to  t,  we  also  have 


(5.12)  «(P,u,9)tl=  ♦  l»tt,0tt.2  ♦  IP„>’  * 


♦  <  w2  *  'f'2  *  "f’t'2 +  *  "tv2** 


4  4 

<  C(»Pfl,u  ,9  I2  +  l  If 1  (0 )  +  sup  (  l  IfV  +  IfV) 

i-0  0<T<t  i-1  z  z  1 


t  A  .  4 

m2 


*  {  ,,°'2  *  J,  ,fi'i  *  ♦  J,  <•>> 


+  CN  (0,t)  . 

By  beginning  with  the  estimates  (5.11)  and  (5.12)  and  combining  the  Lemmas 
5.2  “  5.10  again,  we  can  reach 

(5.13)  sup  («P,u,9l2  +  «u  , 9  I2  +  «u  ,8  I2  +  Ip  I2  + 

0<T<t  q  t  t  2  tt  tt  t  3 


+  ^tt1?  +  /0  ^e,5  +  ,pi;  +  w2*  + 


+  V'Vi ♦  '<£  V2 *  ,(i>tt,2dt  < 
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1 


<  C(IPft,un,9  +  SUp  (lf°l*  ♦  IfJlJ  ♦  I  (If1!! 
0  0  0  4  0<t<t  3  *  1  i=1 


Iffl2)  +  J*  lf°l2  +  Ifjl2  +  If®  I2  +  l  ( 

t  0  i«1 


If1!2  +  If^l2  +  |ftt,-1)dT  +  * 


Therefore,  noting  that 


dp  -  i 

~  -  Pfc  +  P  uJPx  , 


we  have  consequently 


Lemma  5*11.  It  holds  that 


(5.14)  N2(0,t)  <  C(IP0,u0,6o|2  +  sup  (lf°l2  +  lf°«2  + 

+  I  (If1!2  +  If^l2))  +  ^  (lf°l2  +  *f°l2  + 

i=1  0 

,ftt|2  +  ^  (,fi,3  +  ,ft'5  +  ,fJt,-i)dT  +  • 

Now  let  us  establish  the  estimates  for  {f*}. 

Lemma  5.12.  It  holds  that  the  right  hand  side  of 
(5.14)  <  C(IP0,u0,90l2  +  (N(0,t)  +  l*l5)N2(0,t)). 

Here  we  note  that  Lemmas  5.11  and  5.12  easily  imply  the  desired  a  priori 
estimates,  in  fact,  we  may  choose  so  small  that 

N2(0,t)  <  C,P0,u0,9q for  N(0,t),  l*»5  <  ei  . 

Proof  of  Lemma  5.12.  Because  there  are  too  many  terms  to  estimate,  we 
only  give  proofs  for  a  few  examples.  The  remaining  terms  can  be  proved  in  the 
same  way. 

Consider  f®,  for  example,  as  follows: 
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sup  l(p-p-p)u^  -  P  u^ I? 
T  xj  j 


<  sup  (Ip'^ul^  +  Ip-pl^lulJ  +  iDPlJlulJ) 
^  J  4  3  4  3  3 


<  C  E0(N(0,t)  +  )K  (0#t)  , 


/  —  v  «i  v  4  2 

l(p-p-p)uJ  -  P  uJl,  dT 
x .  x .  4 

0  j  } 


<  /  Ipl^lul^  +  Ip-pl^lul^  +  iDpl^lul^  dT 
g  4  3  4  5  4  4 


<  C  E0<N(0,t)  +  »♦!  )N  <0,t>  , 


and  so  on.  Proceeding  in  this  manner  proves  Lemma  5.12.  Finally,  in  order  to 
prove  Theorem  1.1,  it  suffices  to  show  the  exponential  decay  of  the  solution. 
To  do  that,  we  may  estimate  the  equations 


L0(.W«..a))  -  .“V  *  ®atP  , 

u 


_  i .  CM-  n  at  i  at  i 

L  (e  (P,u,8))  =  e  f  +  ae  u  , 


_ 4 ,  at.  at.4  ,  „  at 

L  (e  (u,9))  =  e  f  +  ae  8  , 


in  the  same  way  as  (5.11)  and  taking  a  suitably  small. 
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